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We introduce a family of states, the fPEPS, which describes fermionic systems on lattices in 
arbitrary spatial dimensions. It constitutes the natural extension of another family of states, the 
PEPS, which efficiently approximate ground and thermal states of spin systems with short-range 
interactions. We give an explicit mapping between those families, which allows us to extend previous 
simulation methods to fermionic systems. We also show that fPEPS naturally arise as exact ground 
states of certain fermionic Hamiltonians. We give an example of such a Hamiltonian, exhibiting 
criticality while obeying an area law. 



I. INTRODUCTION 

Understanding the behavior of correlated quantum 
many-body systems is one of the most challenging prob- 
lems in various fields of physics. For spin systems on a 
lattice with local (i.e., short-range) interactions, power- 
ful methods have been developed in recent years. They 
rely on families of states which, on the one hand, de- 
pend on very few parameters and, on the other, ap- 
proximate the quantum state of the spins in thermal 
equilibrium. In one spatial dimension. Matrix Product 
States (MPS) [1^ (which underly [H |3] the successful 
Density Matrix Renormalization Group algorithm [H |^ ) 
provide a good approximation to the ground state of 
any gapped local Hamiltonian. Projected Entangled Pair 
States (PEPS) [6,l^ (cf. also [SIS]), which naturally ex- 
tend MPS to higher spatial dimensions, approximate spin 
states at any finite temperature [TO^, and have been suc- 
cessfully used to simulate spin systems which cannot be 
dealt with otherwise plT[T4] . 

Fermionic quantum many-body systems are central to 
many of the most fascinating effects in condensed matter 
physics. In one spatial dimension, it is possible to adapt 
the methods based on MPS to such systems thanks to 
the Jordan- Wigner transform, which maps fermions into 
spins while keeping the interactions local. In higher di- 
mensions, however, this is no longer possible: fermionic 
operators at different locations anticommute, which ef- 
fectively induces nonlocal effects when mapping fermions 
to spins. Thus, the use of PEPS to describe fermionic 
systems is no longer justified (see however [151 HI] foi" 
different approaches). 

In this article we introduce a new family of states, 
the fermionic Projected Entangled Pair States (fPEPS), 
which naturally extend the PEPS to fermionic systems. 
According to their definition, fPEPS are well suited to 
describe fermionic systems with local interactions. They 
can be, in turn, efficiently described in terms of standard 
PEPS at the prize of having to double the number of pa- 
rameters. This automatically implies that the algorithms 
introduced to simulate ground and thermal states, as well 
as the time evolution of spin systems using PEPS [6 , 7 , 
can be readily adapted to fPEPS. We also show that cer- 
tain fPEPS are exact ground states of local fermionic 



Hamiltonians, in as much the same way as PEPS are for 
spins [17 . In particular, we give the explicit construction 
of a Gaussian Hamiltonian which has a fPEPS as exact 
ground state. Remarkably, the state is critical, i.e. gap- 
less with polynomially decaying correlations, yet obeys 
an entropic area law [^, in contrast to what happens 
with other free fermion systems [19 . 

We have organized this paper as follows. First, we 
will briefly review PEPS and explain why they are well 
suited to describe spin systems with local interactions in 
thermal equilibrium. Then, we will construct the family 
of fPEPS following the same idea. We will then consider a 
subfamily of fPEPS for which we can build local "parent" 
Hamiltonians, i.e. those for which they are exact ground 
states. Finally, we will give a particular example which 
presents criticality. For the sake of simplicity, we will 
concentrate on two spatial dimensions. 

II. CONSTRUCTION OF PEPS 

For simplicity, let us consider a 2D lattice of = 
Nh • Ny spin 1/2 particles, with states |0) and |1). To 
each node of coordinates (/i, v) we associate four auxil- 
iary spins, with states \n) (n = 0, . . . , I) — 1), where D is 
called bond dimension. Each of them is in a maximally 
entangled stateV |n, n) with one of its neighbors, as in- 
dicated in Fig.fTkT The PEPS |^) is obtained by applying 
a linear operator ("projector") to each node that maps 
the auxiliary spins onto the original ones. This operator 
can be parametrized as 

D-l 1 

P(Kv)= E T.^B(!,^-))fluJk){lr,u,d\. (1) 

l,r,u,d=0 k=0 

Let us now explain why PEPS are well suited to de- 
scribe spins in thermal equilibrium in the case of local 
Hamiltonians, H = "^hx. For simplicity, we will as- 
sume that each hx acts on two neighboring spins. We 
first rewrite the (unnormalized) density operator e~^^ = 
^^b[\x) (xI]? where |^) = e~^^/^ 1 \x)ab ^ purifica- 
tion [20^ and \x)ab ^ pairwise maximally entangled state 
of each spin with another one, the latter playing the role 
of an environment. We will show now that |^) can be 
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FIG. 1: a) Construction of a PEPS in two dimensions. The 
balls joined by lines represent pairs of maximally entangled 
D-dimensional auxiliary spins, which are then mapped to the 
physical spins (red), as illustrated by the light blue spheres, 
b) Why PEPS approximate thermal states well: exp[— 
can be implemented using local maps only if an entangled pair 
is available. 

expressed as a PEPS. We consider first the simplest case 
where [hx, hy] = 0, so that |^) = Ux e'^^^/^ ^ 1 \x)ab' 
The action of each of the terms e~^^^/^ on two spins 
in neighboring nodes can be viewed as follows: we first 
include two auxiliary spins, one in each node, in a max- 
imally entangled state, and then we apply a local map 
in each of the nodes which involves the real spin and the 
auxiliary spin, which ends up in |0). By proceeding in 
the same way for each term e~^^^/^, we end up with the 
PEPS description (see Fig.[l]3). This is valid for all values 
of /3, in particular for /3 ^ oo, i.e., for the ground state. 
In case the local Hamiltonians do not commute, a more 
sophisticated proof is required [10]. One can, however, 
understand qualitatively why the construction remains to 
be valid by using a Trotter decomposition to approximate 

e-^^ ^ Yl^^iUx^'^^"^^^^ with M < 1. Again, this 
allows for a direct implementation of each exp[— /3/ia/2M] 
using one entangled bond, yielding M bonds for each ver- 
tex of the lattice. Since, however, the entanglement in- 
duced by each exp[— /3/ia/2M] is very small, each of these 
bonds will only need to be weakly entangled, and the M 
bonds can thus be well approximated by a maximally 
entangled state of low dimension. Note that the spins 
belonging to the purification do not play any special role 
in this construction, and thus we will omit them in the 
following. 



III. CONSTRUCTION OF FPEPS 

We will now extend the above construction to fermionic 
systems, in such a way that the same arguments apply. 
We consider fermions on a lattice, and work in second 
quantization. For a Hamiltonian H = ^hx^ each term 
hx must contain an even number of fermionic operators. 



in order for the Trotter decomposition to be still possible. 
Thus, we just have to find out how to express the action 
of e~^^^ in terms of auxiliary systems. This is very sim- 
ple: one just has to consider that the auxiliary particles 
are fermions themselves, forming maximally entangled 
states, and write a general operator which performs the 
mapping as before. Following this route, we arrive at the 
definition of fPEPS. More specifically, we define at each 
node (/i, v) four auxiliary fermionic modes, with creation 
operators /^(V^), 7(V^)' ^(V^)' respectively. We de- 

fine 

H^n,.) = ;^(l + /5(V«)«Ui,.)) (2) 
V^K.) = ^(l + <5(V„)7?,.„+i)) (3) 

which create maximally entangled states out of the vac- 
uum. We also define the "projectors" 

where a^/j^ is the annihilation operator of the physical 
fermionic mode, and the sum runs for all the indices from 
to 1, with the condition that {u^d^l^r^k) mod 2 = c, 
were c is fixed for each node [26]. The latter is related 
to the parity of the hx and will ensure that the parity of 
the fPEPS is well defined. The fPEPS is then 

1^) = ( n ^{h,v) n ^{v,h)y{v,h))^u^ |vac), (5) 

{h,v) ih,v) 

where the expectation value is taken in the vacuum of 
the auxiliary modes, and |vac) denotes the vacuum of 
the physical fermions. Note that the definition of fPEPS 
straightforwardly extends to systems with both more 
than one physical mode per site and more than one mode 
per bond, as well as to open boundaries or higher spatial 
dimensions. 



IV. RELATION BETWEEN FPEPS AND PEPS 

Next, we will find an efficient description of any fPEPS 
in terms of standard PEPS. With that, one can read- 
ily use the methods introduced for PEPS [6] [7] in order 
to determine physical observables, as well as to perform 
simulations of ground or thermal states, and time evolu- 
tion. We have to identify the Fock space of the fermionic 
modes with the Hilbert space of spins. For that, we sort 
the lattice sites according to M = {v — l)Nh + h and 
associate a\^^ . . . a]y^^|vac) to the spin state . . . , /cat). 
Then we write |^) in that basis, and express it as a PEPS 
in terms of tensors B ([T]). The goal is to find the rela- 
tion between the tensors B (corresponding to the spin 
description) and A (fermionic description). In principle, 
the fPEPS to PEPS transformation can be done straight- 
forwardly by adding extra bonds to the PEPS which take 
care of the signs which arise from reordering the fermionic 
operators; however, this would lead to a linear number 
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of bonds per link and thus to a dimension which is ex- 
ponential in N. Remarkably, it is possible to express 
every fPEPS as a PEPS by introducing only one addi- 
tional bond per horizontal link as follows: Replace each 
fermionic bond by a bond of maximally entangled spins, 
adding one additional horizontal qubit bond everywhere 
except at the boundaries (see Fig. [2|. This means that 
the tensor B will have now two more indices, say I' and 
r', which are associated to those new bonds. Then, we 
find the relation 



(6) 



for = 1, while for h > 1 we have 

{^h,v)\Jrr'ud ~ C^) 
(-l)^<''-'^''"'''''''^(A.„)lL(-l)*''<5;',(r'+«+d)mod2 , 

where f(^h,v) {k^ r) is a function which only depends 

on the local indices, and = for = Nh- 

Let us briefly explain how to obtain this result. Con- 
sider an fPEPS of the form ([5| which we want to bring 
into the normal ordered form by commuting the fermionic 
operators. To this end, we perform the following three 
steps on the total projector ]\Q{h,v)i observing that 
local sign contributions can be absorbed in the ten- 
sors {A{h,v))f}ud' First, commute all physical modes to 
the left. This results in a factor (— where 
P — v) ^ih,v) is the parity of the fPEPS; since the lat- 
ter is fixed, this yields a global phase. Next, contract the 
horizontal bonds: The non-boundary bonds only yield 
local contributions, while the horizontal boundary bond 
on any line v gives a contribution (—1)^(1,^)^(1,^) ^[^i^ 
n(/i, v) = '^jyh{u{j,v) '^^ij^v))' Finally, contract the ver- 
tical bonds, proceeding columnwise from h = 1. For each 
bond between {h^v) and {h,v + 1) this gives a sign con- 
tribution ^-i)d^h,.)Tiih,v). ^^^^ parity of the 

bonds this holds even for the bonds across the boundary. 
Thus, all signs can be computed if the respective par- 
ity Il{h,v) is available at each site, which is achieved by 
the additional bonds passing this information to the left. 
Note that the same proof applies to open boundaries, as 
well as systems with more physical or virtual modes per 
site, without the need for further extra bonds to compute 
Il{h,v). Similarly, one can derive a corresponding result 
for higher dimensions. 




FIG. 2: Every fPEPS can be represented as a PEPS at an 
extra cost of at most one additional bond per link (shown for 
a 3 X 3 PBC lattice). 



thus completely characterized by their covariance ma- 



trix 

.(2) 



kl - where c\'^ = a] + and 

c- ' = {—i){al — ai) are Majorana operators. We will 
now introduce Gaussian fPEPS, which we then use to 
show that fPEPS naturally appear as ground states of 
free local Hamiltonians. The techniques used here follow 
closely the corresponding methods for bosons introduced 
in [21 . 

Gaussian fPEPS are obtained by restricting the map 
Q to be Gaussian {H and V are already of that form). 
Those transform Gaussian states into Gaussian states, so 
that they can be characterized through the map Fin 
Fout- The most general (pure) map can be written as [22] 
Fout = B{D - Fin)-^5^ + A with 



.(1) 



G 



A B 
-B^ D 



-G^ , GG^ 



(8) 



We denote the CM of the translationally invariant states 
of the virtual modes by Fin = ®^h,v where ujh,v is the 
CM of the maximally entangled horizontal resp. verti- 
cal bonds. Then the desired family of states can be ob- 
tained by applying the same Gaussian map to each node 
ft = {h^v) of the lattice: Q = 0fj^, where QQ^ = 1. 
Due to translational invariance, Fout can be conveniently 
expressed in Fourier space, Fout = ®0rout(</^)7 with 

fout = B[D - co{$)]-'BT + A, where = ^) 
is the reciprocal lattice vector. As we show in Appenclix 
[a1 it is straightforward to see that the ^-dependence of 
cj(^) yields 



V. FERMIONIC GAUSSIAN STATES AND 
PARENT HAMILTONIANS 

Fermionic Gaussian states [22 (also known as quasi- 
free states) constitute an important subclass of states, as 
they appear as ground and thermal states of quadratic 
Hamiltonians, corresponding to free fermion or BCS 
states. These states can be written as an exponential 
of a quadratic form in the fermionic operators, and are 



rout((/>) 



(9) 



1 



d((/)) 



Re{qm -lm{q{^)) p{c^) \ 

lra{qm -p{$) Re(g(^)) 

-p{$) -lm{qm -Re{qm / 



with p, q and d low-degree polynomials in (/>; in partic- 
ular, d{(j)) = det[D — oj{(p)]. Now define the Hamiltonian 
H = i hkiCkCi^ where h is defined through its Fourier 
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transform h{(j)) = (i(^)f out(^)- H has Pout as its ground 
state, since rout(</>) and h{(j)) are diagonal in the same 
basis, and unless H is gapless — corresponding to zeros of 
d{(j)) — the ground state is unique. Moreover, since the 
degree of p and q is bounded by twice the number of 
virtual modes per site, it follows that H is local. 

Example Let us now give an example of a local Hamil- 
tonian which has an fPEPS as its exact ground state. We 
present only the main results in the main text, and re- 
fer the reader interested in the details to Appendix [Bj 
We choose the (translational invariant) fPEPS projector 

Q _ g(iQ!+/3)(-7+i5)+Q;/3+75+a'''(-iQ!-/3-7+i5) ^J^j^h yiclds 

p{(j)) / d{(j)) = (sin^i — sin^2)/(— 1 + sin0isin02) and 
q{<j))/d{(p) = cos cos (/)2/(— 1 + sin sin (/)2). The re- 
sulting parent Hamiltonian is (Fig. Ist 



iJcrit — '^'^ ^lh,v)^lh,v^l) ^lh,v)^lh^l,v) 

and Ny odd, which will ensure that the ground state 
is unique. By Fourier transforming To^t{(j)) into po- 
sition space, one obtains that the correlations of Ma- 
jorana operators of equal (different) type at distance 
(ni, 712) scale asymptotically as the real (imaginary) part 
of K{ni,n2) = (ni +3 + m2)/(ni + 1 + m2)^ for ni +n2 
odd (even) and vanish otherwise (Fig. [sJd). Notably, the 
ground state possesses correlations that decay as power 
laws and the Hamiltonian is gapless in the limit N ^ 00. 
In fact, our example provides us with a critical fermionic 
system obeying the area law, which directly follows from 
the fact that its ground state is a PEPS with bounded 
bond dimensions. Note that, although i^crit is not par- 
ticle conserving, it can be converted into a particle con- 
serving one via a simple particle-hole transformation in 
the B sublattice. This new Hamiltonian possesses a spec- 
trum with a Dirac point separating the modes with pos- 
itive and negative energies. Thus, the Fermi surface has 
zero dimension, which explains why our results do not 
contradict the violation of the area law expected for free 
Fermionic systems [19]. 

In summary, in this work we have introduced fermionic 
PEPS (fPEPS) which are obtained by applying fermionic 
linear maps to maximally entangled fermionic states 
placed between nearest neighbors. This construction re- 
sembles the construction of PEPS and is well suited to 
describe ground and thermal states of local fermionic 
Hamiltonians (both free and interacting), in the same 
way as PEPS are suited to describe ground states of lo- 
cal spin systems. We have then shown how fPEPS can be 
transformed into PEPS at the cost of only one additional 
bond providing an explicit mapping for the correspond- 
ing tensors. This also demonstrates the use of fPEPS for 
numerical simulations. Further, we have investigated the 
role of fPEPS as ground states of local Hamiltonians. To 
this end, we have introduced Gaussian fPEPS and shown 
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FIG. 3: a) Hopping terms in the Hamiltonian. b) Exact 
value and asymptotic scaling for the correlations in direction 
of the axis (red), along the diagonal (green) and along the 
direction (n, 2n) (blue) [cf. c)], for Majorana operators of the 
same type (top) and different type (bottom). 



that they naturally arise as ground states of quasi-free 
local Hamiltonians. Finally, we have used these tools to 
demonstrate the existence of local free fermionic Hamil- 
tonians which are critial without violating the area law. 
Note added: After submission of this manuscript several 
algorithms based on fPEPS have been developed and ap- 
plied to interacting fermions [23H25] . 



Acknowledgments 

This work has been supported by the EU projects 
QUEVADIS and SCALA, the DFG Forschergruppe 
635, the excellence cluster Munich Advanced Photonics 
(MAP), SFB Foqus and the Elite Network of Bavaria 
programme QGGG. 



Appendix A: Gaussian fPEPS and Parent 
Hamiltonians 



In this Section we present the details that lead to 
Eq. (9). Recall, that we apply the translationally in- 
variant channel, Eq.(8), to the translationally invariant 
input state Fin = ^0Jh,v The structure of the prob- 
lem suggests an approach in Fourier space, and we in- 
troduce the Fourier transform of the the mode operators 



'^^'^ ffi-, where / is either a physical 

or virtual mode, and (j) = (^f^, ^jf^) is the reciprocal 
lattice vector. Now we consider the CM of the output 
state in the qp-ordered form, i.e. we write 



Tout — 



out out 
T.(2,l) T-(2.2) 



(Al) 



where F, 



(r,s) 



-out - r,s = 0,1. The translation- 

ally invariant construction is reflected in the fact that the 
blocks Fq^'^^ are circulant matrices. Hence, they all can 
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be diagonalized simultaneously by a Fourier transforma- 
tion T . The Fourier transform of Fout, Gout = -^rout-^^ 
has diagonal blocks 

= ■^(^[c('-\cW]>-Ft = (j[d»,d»t]) 



diag W), 



where ^(''•^H'/') ^ C are the eigenvalues of the blocks 



^iut^- The operators d^^^ are the Fourier transformed 



7(r) 



Majorana operators, (i^^ = Xln 
while the Majorana operators in the reciprocal lattice 
space are given by c>^ = ^"x + ^d>^ = (~^)(^^: ~ ^c?)' 



with CM {ttt^)$^^$^ = Both represen- 

tations are linked via a unitary transformation. In the 
following we make use of Gout to derive properties of 
Tout- To this end, we regroup the modes such that 
Gout = ®fGoutW is a direct sum of blocks correspond- 
ing to the same lattice vector, i.e. we write 



^(2,1) ^(2,2) (^) 



(A2) 



Since Fout is antisymmetric and corresponds to a pure 



state, i.e. FH„ 



-1, and the Fourier transformation is 



unitary, we find that Gout(0) can be written as 



Gout((/>) 



1 ( ipicl^) 



(A3) 



where <7((/>), <i((/>) G R. To obtain more information 
on theses functions, we use the fact that the channel £ 
describes a translationally invariant map. This implies 
that the blocks A, B and D are block diagonal, and thus 
commute with the Fourier transform. Hence, 

Gout = ^rout^^ = B{D - G-M)~^B^ + (A4) 

where Gin = J^G-.r^P . We use that [D - Gin(^))"^ = 
adj(I) — Gin(^))/ dei{D — Gin(0)) where adj denotes the 
adjugate matrix, and we define d{(j)) = dei{D — Gin (</>))• 
As Fin is the covariance matrix of a system of maximally 
entangled states between nearest neighbors, its Fourier 
transform Gin((^) is built out of terms of the form e*'^^ ^ 
only. Thus, d{$) = det(D - Gin(^)) and adj(D - Gin(^)) 
are polynomials of low order in (/)i,2- As B and A are 
local operators, we see that p and q are polynomials of 
low degree as well. These results lead to the Fout given 
inEq. (9). 



Appendix B: Example of a critical fPEPS 

Like every Gaussian map the projector Q = 

g(ia+/3)(-7+i5)+a/3+75+at(-ia-/3-7+i<5) described as 



a channel of the form given in Eq. (8), where 



B = 



D = 
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and A = 0. Using this representation, a straightfor- 
ward calculation shows that the functions p, and 
d defined in Eq. (9) are of the form / d{(j)) = 
(sin(/)i — sin^2)/(— 1 + sin (/)i sin (/)2) and q{(j)) / d{(j)) = 
cos cos (/)2/(— 1 + sin sin ^2)- 

Note that the success probability of the PEPS projec- 
tion is related to the absolute value of [22] 



det{D — Gin(0)) cx (—1 — sin 01 sin 02 ) sin 



2 sin^ f 



which means that the fPEPS has zero norm (i.e., is not 
properly defined) for = (|, ^) and = (^, |), as 
well as for 01 = 0, tt and 02 = 0, tt. The former condition 
implies that the state is not defined if the lattice size is a 
multiple of four in both directions. This condition can- 
not be removed, since it is inherent to the way the critical 
model is constructed - these are exactly the zeros of d{(j)). 
The other zeros, however, cancel out in ( A4), and it turns 



out that one can modify the fPEPS construction to have 
nonzero norm in those cases, without changing the CM 
of the state itself. This can be seen by expressing the 
virtual fermions in terms of two Majorana modes: one 
finds that only one of these modes per virtual fermion is 
connected to the physical fermion by the PEPS projec- 
tor, while the other is only perfectly correlated with the 
corresponding Majorana mode of the opposite fermion. 
Thus, these "unused" Majorana modes from perfectly 
correlated loops around the torus, which make the state 
vanish for even loop sizes due to the fermionic statistics. 
By properly modifying the bond across the boundary in 
the unused Majorana mode one can prevent the state 
from vanishing without affecting the fPEPS itself, which 
is still described by (A4). 



Let us now show that the system is critical by deriving 
the asymptotic behavior of the correlation functions p 
and q in position space. For large systems we can replace 
the discrete Fourier transform by a continuous one. Let 
i = p^q and define 



^ni,n2 — (^^(l^l)^(niln2)) 



1 



(2^)2 70 6^(01,02) 

where = 1 and yq = 2 for p- and g-correlations re- 
spectively. We make the substitution z = e^^'^ so that 
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dz = izd(j)i^ and arrive at 

p{z,(t)2) _ . Z 



1 — 2iz sin 02 



d{z^(j)2) 2iz — {z'^ — 1) sm.(j)2^ 
q{z,(j)2) . (2;2 + l)cos02 



(Bl) 



(B2) 



d{z^(j)2) 2iz — {z'^ — 1) sm(j)2 
Then ^ni,n2 can be written as 

Cn„„. = ^ r dcP2e'--^^li^l„^{cP2), (B3) 

where C is the closed loop on the unit circle in the 
complex plane. Since d{z^(j)2)~^ has poles at z± = 
i (1 ± I cos 02 1) / sin 02 and p as well as q are holomor- 
phic at z±^ the integral Inx,n2{^2) is proportional to the 
residue within C according to the residue theorem. As 
only Z- lies within the unit circle we obtain 



(B4) 



Next, we calculate the p-correlations. With the help of 
(Bl ) we obtain 



lip) u,) = J- / dzi 1 

= i-+i(i-|cos</>2ir^^^ 



ni + 1 • 



(sin 02) 

From the symmetry 4^^(02 +^) = (-Ij^^+'^Sn^ (02) 
and the fact that for 02 G [— 7r/2,7r/2] we have | cos 02 1 = 
cos 02, implying 



-+1^5^/ tan 



1 — COS 02 



ni + 1 



(B5) 



we can conclude that 



Pni ,n2 
1 

2^ 



(l-(-l) 



ni+n2 



)2Re 



Jo 



To obtain that result, we have further made use of the 
relation /i^n2(-02) = (-Ij^'^+^^^n^ (02), so that for m 
even (odd) /n^,n2(02) is an odd (even) function, and only 
the sine (cosine)-part of the exponential e*"^2 02 gives a 
non- vanishing contribution. Following a similar strategy, 
one can derive 



To prove criticality we are interested in the asymptotic 
behavior of the integral 

Jo 

The correlations are symmetric under the exchange of 
ni and n2. This follows from translational invariance 
and can also be seen directly from the form of p(0i, 02) 
and g(0i,02). Hence, to determine the asymptotic be- 
havior, we can assume wlog. ni ^ 1. In this limit, 
the absolute value of Ini,n2{4'2) attains its maximum for 
02 = ±arccos(l/n — 1) ^ 7r/2. We rewrite 

where the functions c and t are given by c(02) = i^cos^ 



and t(02) = log (^tan^y Next, we expand c(02) and 
t(02) around 7r/2: = _ (.^^ - f ) + {<j^2 - f )' + 

O ((02 - f )^), log (tan f ) = - I) + i (02 - 1)^ + 
O ((02 - f )^). Substituting 02 ^ 02 - f the integral at- 
tains the form Jni,n2 = ^^^^^^^"^ /^7,/2 ^^2-^(^1, ^2, 02) 
with kernel 

J(ni,n2,02) 

= e^"^^2 (-02 + 4) e(-i+i)^2+^2/6(ni+i)(i ^ 0{4)). 



We use Jni,n 
-7r/2 



/_^d02J(ni,n2,02) 



<^02^(^i7 ^2, 02)5 and obtain 



3 + ni + in2 



1 1 

(l + ni+m2)^ ' ^ V^i'^2 



#2^(^1,712,02) = , ^ + C> ( — , — I , 



while the second integral can be bounded by 



-7r/2 



(i02J(ni, 712,02) 



-7r/2 



This gives rise to only an exponentially small correction 
that can be neglected in the asymptotic limit. Summariz- 
ing, we see that the p-correlations are non- vanishing only 
for ni + 7i2 odd, while g-correlations are non-vanishing 
only for rii + n2 even: 



Qni ,ri 



(l + (-l)"i+"=^)2Re 

27r 



Jo 



Pn,,n, - (l-(-ir+"^)Re(^^ 

Qni ,n2 



3 + ni + in2 



(1+711+ ^712)^ 



(l + (-ir+-)Im(-^±^^l±^). 
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